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TOWARDS A MORE GENERAL TYPE OF UNIVARIATE CONSTRAINED 
INTERPOLATION WITH FRACTAL SPLINES 

A.K.B. CHAND, P. VISWANATHAN, AND K.M. REDDY 


Abstract. Recently, in [Electronic Transaction on Numerical Analysis, 41 (2014), pp. 420- 
442] authors introduced a new class of rational cubic fractal interpolation functions with linear 
denominators via fractal perturbation of traditional nonrecursive rational cubic splines and 
investigated their basic shape preserving properties. The main goal of the current article is 
to embark on univariate constrained fractal interpolation that is more general than what was 
considered so far. To this end, we propose some strategies for selecting the parameters of the 
rational fractal spline so that the interpolating curves lie strictly above or below a prescribed 
linear or a quadratic spline function. Approximation property of the proposed rational cubic 
fractal spine is broached by using the Peano kernel theorem as an interlude. The paper also 
provides an illustration of background theory, veined by examples. 

Keywords. Rational fractal interpolation function; Constrained interpolation; Peano Ker¬ 
nel; Convergence. 

AMS subject Classifications. 28A80; 65D05; 65D07; 41A29; 41A30; 41A25. 

1. Introductory Remarks 

Interpolation, which deals with the construction of a function in continuum from its avail¬ 
ability in a finite set of points, is a rather old problem dating back to ancient Babylon and 
Greece. In view of its increasing relevance in this age of ever-increasing digitization, it is quite 
natural that the subject of interpolation is receiving more and more attention. Consequently, a 
multitude of different interpolation schemes are being developed and reported in various places 
in the literature. However, in times where all efforts are directed towards construction of smooth 
functions, the fact that many experimental and natural signals are “rough” having a dense set 
of nondifferentiable points or even nowhere differentiable was easily ignored. 

In 1986, Barnsley ||Tj introduced the idea of fractal interpolation function (FIF) aiming at 
data that require a continuous representation with “high irregularity”. In words of Barnsley: 

.(these functions) appear ideally suited for approximation of naturally occurring functions 

which display some kind of geometrical self-similarity under magnification, for instance, profiles 
of mountain ranges, tops of clouds and horizons over forests, temperatures in flames as a function 
of time, electroencephalograph pen traces and the minute by minute stock market index. In 1989, 
Barnsley and Harrington [2] provided a seminal result on differentiability of fractal functions, 
one which initiated a striking relationship between fractal functions and traditional nonrecursive 
interpolants, and lead to the creation of a subject which is now called fractal splines. Since then, 
many researchers have contributed to the theory of fractal functions by constructing various 
types of FIFs, including Hermite or spline FIFs (see, for instance, 0HS]), hidden variable FIFs 
(see, e.g., [31 0 0), multivariate FIFs (see, e.g., Ei ini)- An operator theoretic formalism for 
fractal functions which enabled them to find extensive applications in various classical branches 
of mathematics was introduced and popularized by Navascues mm- 

In spite of about a quarter century of its first pronouncement, FIFs were extensively investi¬ 
gated only in terms of their algebraic and analytical properties such as fractal dimension, Holder 
exponent, differentiability, integrability, stability, and perturbation errors. Thus, these studies 
completely ignore one of the important properties of an interpolation scheme, namely, preserving 
shape inherent in the data. Constrained control of interpolating curve is a fundamental task 
and is an important subject that we face in applications including computer aided geometric 
design, data visualization, image analysis, cartography. However, fractal functions are not well 
explored in the field of constrained interpolation. Motivated by theoretical and practical needs, 
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the authors have initiated the study of shape preserving interpolation and approximation using 
fractal functions; see, for instance, [T9l H8] . However, these researches concern about preserving 
the three basic shape properties, namely, positivity, monotonicity, and convexity. On the other 
hand, there are practical situations wherein interpolating curves that lie completely above or 
below a prefixed curve, for instance, a polygonal (piecewise linear function) or a quadratic spline 
are sought-after. The current article can be viewed as a contribution in this vein. The impetus 
for the research direction reported here grew out of the works of Qi Duan et al. mm- 

The structure of this article is as follows. We assemble some relevant definitions and facts 
concerning fractal interpolation in Section [2] In Section [3] we briefly recall the rational cubic 
fractal spline studied in [18]. Section [4] is devoted to establish convergence of the rational cubic 
fractal spline. Strategies to select the free parameters involved in the rational cubic FIFs so that 
they render interpolating curves that lie below or above a prescribed spline (piecewise defined 
function) are enunciated in Section [5] In section [(3 we provide some numerical examples. 

2. Univariate Fractal Function: An Overview 

This section targets to a constellation of a few rudiments of fractal interpolation theory. For 
a detailed exposition, the reader may consult the references mm- 

The following notation and terminologies will be used throughout the article. The set of real 
numbers will be denoted by R, whilst set of natural numbers by N. For a fixed N £ N, we shall 
write Nat for the set of first N natural numbers. Given real numbers x\ and xat with x± < xjv, 
we define C k [xi, Xat] to be the space of all real-valued functions on [xi,Xjv] that are fc-times 
differentiable with continuous fc-th derivative. 

Let N > 2 and {(xj, yi) £ R 2 : i £ Nat} denote the Cartesian coordinates of a finite set of points 
in the Euclidean plane with strictly increasing abscissae. Set I = [xi,Xat] and A = [xj,Xj+i] 
for i £ Nat-i- In fractal interpolation, a continuous function g : I —y R satisfying g{xi) = yi 
for all* £ Nat and whose graph G(g) is a fractal (self-referential set) in the sense that G{g) is 
a union of transformed copies of itself is sought for. There are two approaches for constructing 
fractal interpolation functions. First method introduced by Barnsley characterizes the graph of 
FIF as an attractor of specifically chosen iterated function system, hence the graph of the fractal 
interpolant can be approximated by the “chaos game” algorithm [¥j. In what follows, we supply 
the second approach, to wit, the construction of a fractal function as the unique fixed point 
of Read-Bajraktarevic operator defined on a suitable function space, which was popularized by 
Massopust [T2l . 

Suppose Li : I —»• 7j, i £ Nat-i be affine maps satisfying 

Li{x l) — Xi, Li (xat ) — Xi- (-1- 

For i £ Nat-i, let 0 < ki < 1 and Fi : I x R — > R be functions that are continuous in first 
argument and fulfill 

Fi{xi,yi) =Vi, F^xniVn) = Vi+i, \Fi(x, y) — F i (x 1 y*)\ < k t \y — y*\. 

Define 

C*(I) = {he C(I) : h(x i) = yi,h(x N ) = y N }, C**{I) = {h £ C(I ) : h{xi) = y^i £ Nat}. 

It is valuable to note that both C*(I) and C**{I ) are closed (metric) subspaces of the Banach 
space (C(J), ||.||oo). Define T : C*(I) —>■ C**{I) C C*{I) via 

Th{x) = F i (L~ 1 {x), h o L“ 1 (x)), x £ h, i £ Njv-i- 

The mapping T is a contraction with contraction factor k := maxjfci : i £ Nat-i}, and the fixed 
point g of T interpolates the data {(x*, y*) : i £ Nat}. Consequently, g enjoys the functional 
equation 

g{Li{x)) = Fi(x,g(x )), x £ /, i £ Nat-i. 

By the Banach fixed point theorem, it follows that g can be evaluated by using 

9 = lim T k (g 0 ), 

k—too 
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where T k denotes the fc-fold composition of T, and go G C*(I) is arbitrary. Further, by the 
collage theorem, for any h G C*(I) 

||/i-0||oo< J^-Wh-Thiu 
1 — c 

If Wi : I x R —> Ii x M is defined by Wi{x,y) = (-Lj(x), Fj(x, y)) V i G Njy-i, then 

G{g) = Ui & n N _ 1 w i (G(g)). 

Next let us point out the well-known fact that the notion of fractal interpolation can be applied 
to associate a family of fractal functions with a prescribed function / G C(I). This was observed 
originally by Barnsley and explored in detail by Navascues. Consider the maps 

Li(x) = aiX + Cj, Fi(x,y ) = a z y + f o L z (x) - aiib(x), i G Njv-i, 

where f b : I — >■ R is a continuous function interpolating / at the extremes of the interval 
and cti are real parameters satisfying |ai| < 1, termed scaling factors. The corresponding FIF 
denoted by f& b = /“ is referred to as a-fractal function for / (fractal perturbation of /) with 
scale vector a, base function b and partition A = {xi, x 2 , ■ ■ ■ , Xn} of I with increasing abscissae. 
Note that here the interpolation points are {(xj,/(xj)) : i G Nat}. The function /“ may 
be nondifferentiable and its fractal dimension (Minkowski dimension or Hausdorff dimension) 
depends on the parameter a G (—1, l)^' 1 . Further, /“ satisfies the functional equation 

(2.1) f a (x) = f(x) +ai(f a - x G I it i G Njv-i- 

With the notation |a|oo := maxllctjl : i G Nat_i}, the following inequality points towards the 
approximation of / with its fractal perturbation /“. 

(2-2) ||/“-/|[oo< r !£ H-|| / -6|| 0 o. 

1 - |a;|oo 


3. Rational Cubic Spline FIF Revisited 


Let {( Xi,yi ) : i G Nat} be a prescribed set of interpolation data with strictly increasing 
abscissae. Let dt be the derivative value at the knot point Xi which is given or estimated with 
some standard procedures. For i G Njv-i, denote by hi , the local mesh spacing x i+ i — Xj, and 
let ?'i and ti be positive parameters referred to as shape parameters for obvious reasons. A 
CLcontinuous rational cubic spline with linear denominator was introduced in m as follows. 


f(Li(x)) 


(1 - 0) 3 riyi + 0(1 - 9) 2 Vj + d 2 (l - 9)Wj + 9 3 Uy i+ 1 
(1 - 9)n + OU 


i G Nat-i, 


where 

X — X\ 

Vi = (2 n + ti)yi + Tihidi , Wi = (r< + 2U)y i+ i - Uhidi+i, 6 = —--—. 

xn — 

In reference m, the authors observe that for the construction of fractal perturbation f a of the 
rational spline / involving shape parameters, it is more advantageous to work with a family of 
base functions (6j : i G Nat_i} instead of a single base function used in the traditional setting 
(see Eq. (12.111 L Furthermore, for the fractal perturbation /“ to be C 1 - continuous, it suffices 
to choose the scaling parameters ai such that |cq| < cq and base functions {6j : i G Nat-i} so 
that each bi agrees with / at the extremes of the interpolation interval up to the first derivative. 
Among various possibilities, the following choice is motivated by the simplicity it offers for the 
final expression of the desired fractal analogue of /. 


(3.1) 


bi(x) 


Bu( 1 - Of + B 2l 6{ 1 - 9 ) 2 + R 3 ^ 2 (l - 0) + B 4t 6 3 

(1 - 9)n + 9ti 


where the coefficients Bu, - 821 , - 631 , and Ba are prescribed as 

Bu = nyi, B 2i = (2 n + U)yi + rjdi(xjv - aq), 
B 3 i = fa + 2U)y N - Ud N (x N - xi), B 4i = tiy N . 
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Therefore, in view of m » the desired C 1 -continuous rational cubic fractal spline is given by the 
expression 


/“(£<(*)) 


P i (x) = P*(0) 


M(z) = Q*(9) 


OLif a (.x) + where 

Qi{x) 

(yi - onyi)ri{ 1 - 9) 3 + (y i+1 - aiy N )U9 3 + {(2 u + U)yi 

+r i h i d l - ai[(2 n +ti)yi + t^at - Xi)di]}d(l - 0) 2 + {(r-* + 2ti)y i+ i 

-tihidi+i - ai[(ri + 2 U)y N - U{x N - xi)d N ]}9 2 (l - 9), 

(1 - 9)n + 9ti , i G Njv-i, 9 = ———. 

X AT — Xi 


Remark 3.1. A rational cubic spline with linear denominator which is based only on the function 
values is introduced in m, pointing out that in some manufacturing processes the derivatives 
are difficult to obtain. However, as far as we know, it differs from d3j) and its construction only 
in the following way. Consider the given set of data points Ai = { (x. n yi) : i G Nat+i} and 
the subset A 2 := {(xi,?/i) : i G Nat} of interpolation points. Treating di to be equal to the 
chord slope A, = y '+^~ Vt for % g Nat, the rational cubic spline involving only the function values 
discussed in detail in lOi can be deduced at once from (0. Consequently, on similar lines, the 
following expression for the fractal rational cubic splines involving only function values can be 
obtained. 

f a (Li(x )) = aif a {x) + 

Qi{x) 

Pi(x) = (yi - aiyi)ri(l - 9) 3 + (y i+ i - ctiy N )ti9 3 + 

{(n + ti)yt + rryi +1 - aj[(2r, + U)y\ + n(xN - xi)Ai]}0(l - 9) 2 + 

{ (2 ti V z 1 h.jtj Aj_)_r ctj [(2t ? ; -f- Ti'jyN ti (xn *^l)AA/']}d iJ (l 9 ), 

(3.3) Qi(x) = ri(l - 9) +ti9, i G Njv-i, 9 =———. 

xn — xi 

Hence, in principle, the analysis on the fractal spline structure (13.2[) set about to do in the 
subsequent sections can be easily modified and adapted to the rational cubic spline FIF given 
inEq. ([3D. 


4. Convergence Analysis 

The rate at which an interpolant g approaches the unknown function generating the data 
{(xi,yi) : i G N/v}, is perhaps the most important factor deciding the efficacy of an interpolation 
scheme. In [18], it has been established that the rational cubic fractal spline /“ converges to 
$ G C 2 (I) with respect to the C 2 -norm. Since order of continuity of f a is C 1 , it is natural to ask 
whether /“ possesses uniform convergence as the norm of the partition approaches zero, if the 
data generating function is assumed to have a reduced order of continuity, namely C 1 . In this 
section, we answer this in the affirmative. 

Note that the implicit and recursive nature inherent in the description of the fractal spline /“ 
make it impossible or at least difficult to apply the standard tools for error analysis available 
in the classical interpolation theory for establishing its convergence. On the other hand, we can 
base our analysis on the trustworthy triangle inequality 

(4.1) ii$-nioc<ii$-/iioo + ii/-nioc. 

Following (12.21) . it can be deduced at once that 

( 4 . 2 ) II/- nioo < (II/II00+ max INU). 

1 — a|oo *€N n _i 

It can be read from m that 

ll/lloo < \y\oo + ^\d\oo, 

INloo < max{|j/i|, \y N \} + ^j-max{|di|, |djv|}, 
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where |y|oo := {\yi\ ■ i £ Nw}, |d|oo := {|d*| : i £ N N }, h := {hi : i £ Nat_i}, and \I\ := x N - x\. 
Observing that \at\ < at = jh for all i £ Nat- 1 , we obtain 

|o)oo ^ h 
1 oloo |-f| h 

Hence from (14.21) . the perturbation error obeys 

n/-riioo = <m as /i —)> o. 

As an interlude to our analysis, we have the following theorem which deals with the first summand 
of the inequality (14. Ill , namely, the interpolation error of the traditional rational cubic spline / 
(cf. Eq. ©) with respect to the uniform norm. 

Theorem 4.1. Let {( Xi,yi ) : i £ Nat} be a prescribed set of interpolation data generated by a 
function <f> £ C 1 (/) and f be the corresponding traditional nonrecursive rational cubic spline with 
linear denominator. Assume that the derivatives at the knots are given or estimated by some 
linear approximation methods (for instance, arithmetic mean method). Then the local error of 
the interpolation is given by 

|<f>0) - f(x )| < hi CiU^'Hoo, for x £ Ii, 


where for the local variable ip := x h Xi ■ the constant Ci is given by 


Ci = max 
O<0<1 


rjip( 1 - y?) 2 (l + 2 ip) + tjtp 2 { 1 - y>)(3 - 2 ip) 

n{ 1 - ip) +Uip 


Proof. Consider the error function E($;x) = 4>(a;) — f(x) as a linear functional which operates 
on 4). Direct calculations confirm that the linear functional annihilates polynomials of degree 
strictly less than one. By the Peano kernel theorem m- 


(4.3) 


L[<&] = E(<S>\x) = $(x) — f(x) = 


rx i+ i 

L 


&(t)L x [(x - t)° + ] dr , 


where the kernel function is given by 


L x [(x-t)° + ] 


r(r, x) if Xi < t < x, 
s(t,x ) if x < t < Xi+\. 


Here the notation L x is used to emphasize that the functional L is applied to the truncated 
power function 


(x-tT + :={ (*--)" 
v 1 0 if t > x, 


considered as a function of x. Denoting ip := x , Xi , a rigorous calculation yields 

1 - <p)(l - V 2 ) + ti<p( 1 - ip) 2 


r(r, x) = 


(1 - ip)ri + <pti 

rjip 2 (l -ip)+ tjip 2 (2 - ip) 
(1 - ip)n + ipti 


s{t,x) 
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Using the above expression for the kernel function in (14.311 . we obtain the following chain of 
relations: 


|$(a:)-C(x)| = 


< 


< 



l$'l 


r(r, x ) dr — 


Xi+l 


<(t,x) dr 


14)' 


/M f rj{ 1 - y>)(l - p 2 ) + tjip( 1 - pY 
(1 - ip)ri + <pU 
riV? 2 (l - p) + Uip 2 { 2 - p) 


{X - Xi) 


(1 - p)Ti + pti 

W\\h iCi , 


(xi +1 - a:)}, 


where ||.|| denotes the uniform norm on the subinterval Ii = [xi,Xi+\\. This offers the promised 
result. □ 


Moving now to the crux of this section, we have the following theorem, whose proof is imme¬ 
diate from the foregoing discussion and the theorem. 


Theorem 4.2. Let {(a h,yi) : i € Nat} be a prescribed set of interpolation data generated by a 
function $ € C 1 (J). Assume further that the derivatives at the knots are obtained by some linear 
approximation methods. Let /“ be the corresponding rational cubic spline FIF. Then 


I *-/ c 


OO _ 


< 


1 |^|oo 

+c/i||<f ,/ ||oo 


Moo + max{|yi|, |j/jv|} + + \I\ max{|di|, |djv|}) 


where c := max{cj : i £ Nat_i}. 


The following remark highlights the important fact to which the preceding theorem sheds 
light. 

Remark 4.1. If 4> £ C 1 (I) is the data generating function and f a is the rational cubic fractal 
spline corresponding to the data set, then 

\\*-r\\ 0O = 0{h), as /i —» 0. 

Ergo, /“ converges uniformly to the original function, as the norm of the partition tends to zero. 


Remark 4.2. Let us point out, at least for the sake of good bookkeeping, that the application of 
triangle inequality does not imply that the error in approximating a continuously differentiable 
function with the fractal spline /“ is always greater than or equal to the error in approximation 
by its classical counterpart /. We employed the triangle inequality just to demonstrate that the 
fractal spline /“ has the same order of convergence as that of the classical rational spline /. 
Finding appropriate a for which /“ is close to $ is an entirely different problem, which can be 
shown to be a constrained convex optimization problem with the aid of collage theorem [lf)j . 
There is no reason to believe that the solution would always be a = 0, which corresponds to the 
classical interpolant /. 


5. Constrained Interpolation with Rational Cubic Spline FIF 

This section features a systematic discussion of selection of parameters involved in the ratio¬ 
nal cubic spline FIF, focusing on its applicability in a more general constrained interpolation 
problem. To be precise, given a data set {(a;*, yi) : i £ Njv} and a function p (piecewise linear or 
piecewise quadratic with joints at the knots Xi) satisfying yi > p(xi), the problem is to construct 
a rational cubic spline FIF /“ such that f a (x) > p(x) for all x £ I. The reader is bound to 
have noticed that the values of /“ in the subinterval [a;j,Xi+ 1 ] depends on its values on other 
subintervals as well. Hence, obtaining conditions for which /“ lies above a piecewise defined 
function is relatively harder than that in the corresponding classical counterpart. However, this 
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can be remedied by connecting /“ with its classical counterpart f° = /, and performing the 
analysis on / rather than on f a itself. Let us commence by noting that 

f a {x)~p(x) = f a (x) -f(x) + f(x) -p{x) 

> , ^°° M + f(x)-p(x), 

|a|oo - 1 

where 

M = |y|oo + max{|j/i|, |yjv|} + ^(^|d|oo + \I\ max{|di|, \d N \})- 
For brevity, we denote 

AT = | a l°° M. 

I^loo f 

It follows that for f a (x) — p(x) > 0, it suffices to make 

f{x) > p{x) - K. 

Recall that on [xi,Xi + 1 ], f(x) = where with local variable p := 

Ri{x) = (1 - pf'ri'yi + p{l - tp) 2 [{2n +U)yi + Tihidi] + p 2 (l - p)[(n + 2ti)y i+1 

tihidi- fl] T p Ci/z+1 , 

Si(x) = (1 - p)ri + pU. 

Let p be a piecewise linear function with joints at Xi, i G Njv. If Pi = p{xf) and Pi+\ = p(xi+i), 
then we have 

p{x) =pi( 1 - p) +p i+ ip , x G Ii = [xi,x i+ 1 ]. 

Therefore f(x) > p(x) — K for all x G I is satisfied if 

(5.1) Ri(x) - \pi{ 1 - p) + Pi+ip - K]Si(x) > 0, for x G /*, R Njv-i- 

By using the technique of degree elevation, we assert that 

\pi{l - p) +p i+ ip\Si(x) = riPi(l - p) 3 + [r.i(pi +p i+ i) +PiU]p(l - p) 2 

+ [riPi+ i + ti(pi + p i+ i)\p 2 (l -<!>)+ p 3 Up i+ i, 

Si(x) = (l - p)r.i + pti = n{l - p) 3 + (2 n + ti)( l - p) 2 p + (n + 2u)p 2 {i - p) 

+p 3 U. 

In view of the aforementioned pair of equations, the condition m can be recast as a problem 
of positivity of a cubic polynomial, namely, 

n{yi ~Pi + K)(l - p ) 3 + [n(2yi + hidi - p z - p i+1 + 2 K ) + U{yi - pi + A')] ^(1 - p) 2 
+ [ri(y i+ 1 ~Pi+i + K) + ti(2y i+1 - hid i+1 - pi - p i+ 1 + 2K)]p 2 (l - p) 

+ti{y i+ i -p i+ 1 + K)p 3 > 0, for all i G Njv-i- 

Conditions for positivity (nonnegativity) of a cubic polynomial is well-studied in the literature 
(see, for instance, [17j). However, to keep our analysis simple enough, we shall impose conditions 
on parameters so that each coefficient of the cubic polynomial appearing in the above inequality 
is nonnegative. Noting that j/j > pi for alii G N/v and M > 0, we obtain 

Vj - Pj + K > 0 <*=> Moo < ——— 1 for 3=h i (ii) + !■ 

Vj - Pj + M 

The discussion we had until now is epitomized in the following theorem. 

Theorem 5.1. Suppose that a data set {{Xi,yi) : i G N^r}, where yi > pi = p{xi) and p is a 
piecewise linear function with joints at knots Xi is prescribed. Then a sufficient condition for the 
rational cubic spline FIF f a to lie above p is that the parameters satisfy the following inequalities: 

(i) M| < at for i G Njv-ij and |a| OO - min L * : * e Njv }> 

Vi Pi i dVl 

(ii) ri(2yi - p i+ 1 - pi + hfi + 2 AT) + tfyi - Pi + K) > 0, i G Njv-i, 

(Hi) n(y i+ 1 - p i+ 1 + AT) + ti(2y i+1 - p i+1 - Pi - Mi+i + 2A') > 0, i G 
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where 


M = |y|oo + max{|yi|, \y N \} + -(/i|d|oo + |/| max{|di|, |<ijv|}), 


K = 


-M. 


|o)oo 1 

A couple of remarks that supplement the stated theorem are in order. 


Remark 5.1. On similar lines, conditions for /“ to lie below a piecewise linear function p* 
with joints at Xi satisfying p*(xi) < y,; can be derived. By coupling these conditions with 
those prescribed in Theorem 15.11 we obtain strategies for selecting the scaling factors and shape 
parameters so that p*(x) < f a (x ) < p(x) for all x G I. To keep the article at a reasonable 
length, we avoid the computational details here. 

Remark 5.2. Having selected the scaling parameters according to the prescription in Theorem 
IQ the existence of rational cubic fractal spline /“ that lie above a piecewise linear function p 
points naturally to the existence of positive shape parameters r, and ti satisfying the inequalities 
therein. We shall hint on the existence in the following, see also m- Let us rewrite the 
inequalities involving r* and ti (see items (ii)-(iii) in Theorem 15. II) as follows. 

(5.2) A t + A iBi >0, Cj + \Di > 0, 

where A* = |j-, Ai = 2yi - p i+1 - p t + hidi + 2 K, B t = yi - Pi + K, C t = y i+ 1 - p i+ i + K, and 
Di = 2yi + i —pi +1 — pi — hidi+i + 2K. If pi < yi for all i G Nat and the scaling factors are chosen 
such that |a|oo < min { y V Jp.+M : £ Nat} (that is, Bi > 0 , Ci > 0), then a positive A, (that is, 

a set of positive r 7 ; and ti) satisfying the inequalities in (15.21) exists except when A, > 0, Di < 0, 
and AiDi > 6', B t . 


Remark 5.3. If we replace d t with the slopes Aj, then the coefficient of r* in item (ii) of 
Theorem 15.11 reduces to yi — pi + yt +1 — Pi +i + 2K. Hence, in view of item (i), the condition in 
item (ii) is trivially satisfied. Consequently, we obtain the following conditions for the rational 
cubic spline FIF /“ (cf. (13.31) 1 to lie above a piecewise linear function p. 

\ati\ < at, i G Nat-i, |a|oo < min {——— : i G Nat}, 

Vi - Pi + M 

n(yi + i - p i+ 1 + K) + ti{2y i+1 - p i+ i - pi - hi A i+1 + 2 K) > 0, i G Njv-i- 

Next let us consider p to be a piecewise defined quadratic polynomial with joints at Xi, i G Nat. 
If Pi = p(xi), Pi+i = p(x i+ 1 ), and p'{xi) = p[, then we have 

p(x) = pi( 1 - p) 2 + (2 pi +p'ihi)ip( 1 - tp) +p i+ pp 2 1 x G h= [xi,Xi + 1 ]. 

Analysis similar to the case of piecewise linear p establishes the following theorem. 

Theorem 5.2. Suppose that a data set {(Xi,yi) : i G Nat}, where yi > Pi = p{xt) and p is a 
piecewise quadratic function with knots at Xi is given. Then a sufficient condition for the rational 
cubic spline FIF /“ to lie above p is that the parameters satisfy the following inequalities: 

yi — Pi 

(i) |cq| < di fori G Nat-i, and |a|oo < min {- - -A— : i G Nat}, 

Vi - Pi + M 

(ii) ri(2yi - 2pi + hidi - hip\ + 2 K) + U{yi - pi + K) > 0, i G Nat-i, 


(in) ri(y i+1 - p i+1 + K) + ti(2y i+1 - 2p t - Mi+i - hip\ + 2 K) > 0, i G N^-i, 
where 


M 

K 


I U |oo + max{|yi|, \y N \} + -(h|d|oo + |/| max{|di|, |djv|}), 


- 1 


-M. 


Note that if M is very large, then admissible value for |a|oo will be close to zero. Consequently, 
the fractal perturbation /“ will “almost” coincide with /. This is an unpleasant situation as 
far as construction of a C 1 -continuous constrained interpolant with “fractality” in its derivative 







TOWARDS A MORE GENERAL TYPE OF UNIVARIATE CONSTRAINED INTERPOLATION WITH FRACTAL SPLINES 


is concerned, owing to the fact that larger the value of |a|oo (with respect to the interpolation 
step), more pronounced is the irregularity in the derivative of the smooth fractal interpolant. 
Furthermore, in light of Remark 15.21 it can be deduced that there are some cases in which the 
constrained interpolation cannot be solved with the rational cubic spline FIF /“. These bring 
us motivating influence to seek an alternative approach. We shall now provide a quite brief 
deliberation of this alternative approach for the constrained interpolation. The following is the 
key theorem. 

Theorem 5.3. Let f £ C 1 (I) and A := |aq, aq,..., aw} be a partition of I satisfying x± < aq < 
• • • < aw- Let the scaling factors be chosen such that oii £ (0, a,;) and the base functions be chosen 
fulfilling the conditions b^\x i) = f^\xi), &P(aw) = f^\xisr) for j = 0,1 and bi(x) < f(x) 
for all i £ N/v-i, x £ I. Then the corresponding fractal function /“ € C 1 )/) and f a (x) > f(x) 
for all x £ I. 

Proof. Recall that fractal function /“ corresponding to / satisfies the functional equation 
(5.3) f a (Li(x)) = f(Li{x)) + a i {f°‘{x) - bi(x )], x £ I. 

For |aj| < ai and base functions coinciding with / at the extremes of the interval up to the first 
derivative, it follows [18] that /“ £ C(I). 

Further, f a (xi) = f(xi), i.e., (/“ — f)(xi) > 0 for all i £ Nat (in fact, the equality). Note that 

I = U LAI) and /“ is constructed using an iterated scheme according to (15.31) . Whence, to 
iENjv-i 

prove (/“ — f)(x) > 0 for all x £ I it is enough to prove that (/“ — /) > 0 holds good at the 
points on I obtained at (z + l)-th iteration whenever (/“ — /) > 0 is satisfied for the points on 
I at z-th iteration. The aforementioned condition is equivalent to (/“ — f)[Li(x)) > 0 for all 
i £ Nat_i whenever (/“ — f)[x) > 0. We may rewrite (15.31) in the form 

, . (r-f)(L i (x))=a i (r-b i )(x), 

= a z (f a - f){x) + cxiff-bi)(x). 

Assume that the scaling factors are chosen such that cu > 0 for all i £ Njv-i- By our assumption 
(/“ — f)(x) > 0, hence an appeal to (15.41) reveals that (/“ — f)(Li(x)) > 0 is satisfied whenever 
h(x) < f{x) for all x £ I and i £ Nn-i, offering the proof. □ 

The foregoing theorem demonstrates that if a constrained interpolation problem is solved 
with a classical rational interpolant /, then the perturbation process can be so designed that the 
corresponding fractal spline /“ also provides a solution to the same constrained interpolation 
problem. For instance, let / be the rational cubic spline with linear denominator lying above 
a piecewise linear function p. Choose ai £ (0, Oj) and the base functions 6 j = / — 6 ,;, where 
bi are positive functions satisfying U/\xi) = 0, S-^(xat) = 0 for i £ Njv_i, j = 0,1. Since 
positivity preserving Hermite polynomial and rational spline interpolants are well-studied in the 
literature, it is not hard to find a family of functions {&; : i £ Nat_i} satisfying the aforementioned 
conditions. Note that cubic Hermite interpolants may not serve as candidates for bi, as they may 
reduce to zero functions owing to the imposed conditions. The fractal function f a corresponding 
to / with these choices of scaling parameters and base functions will lie above p. Thus, by dint 
of a more careful choice on base functions bi, the conditions on a, can be made relatively weaker, 
to wit, ai £ ( 0 ,aj) for all i £ Njv-i- 


6. Numerical Examples 

In this section, the developed conditions on parameters are leveraged to obtain rational cubic 
spline FIFs that lie above a prescribed piecewise linear function (polygonal) or quadratic func¬ 
tion. Consider the data set {(aq, yi, di),i = 1,2,3,4,5} = {(0,18, —4.02), (3,10, —1.31), (7,12, —0.36), 
(10,9,0.2), (15,20,4.2)}. Note that the prescribed data set lies above the piecewise linear func¬ 
tion p with nodes at {(0,12), (3,4), (7,10), (10,4)} given by 

— |a; + 12 if 0 < x < 3, 

|a; — \ if 3 < x < 7, 

-2a;+ 24 if 7 < x < 10, 
lx - 10 if 10 < x < 15. 

5 — — 
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FIF. 




(b) Effect of change in a in (c) Classical rational cubic 
Fig. 1(a). spline FIF. 


Figure 1. Rational cubic spline FIFs lying above a piecewise linear function p. 


Let the derivatives at knot points be {—4.02, —1.31, —0.36, 0.2,4.2}. Suppose that, due to some 
reasons, perhaps for a valid physical interpretation of the underlying process, a constrained 
interpolant lying above the following piecewise linear function p is required. For brevity, let us 
represent the parameters, namely, the scaling factors on satisfying |a,;| < < 1, and the positive 

shape parameters r,; and t, as vectors denoted by a, r and t in the four dimensional Euclidean 
space R 4 . The details of the scaling vectors and shape parameters used in the construction of 
constrained rational FIFs are provided in Table [1] for a quick reference. 


Table 1. Scaling vectors and shape parameters used in the construction of 
rational cubic FIFs. 


Figs. 

Scaling Vectors a 

Shape Parameter Vectors r, t 

Fig. (Ha) 
Fig. dlb) 
Fig. die) 
Fig. 2(a) 
Fig. 2(b) 
Fig. 2(c) 

(0.010,0.020,0.030,0.333) 
(0.027,0.027,0.027,0.024) 
(0,0,0,0) 

(0.012,0.013,0.040,0.005) 

(0.040,0.039,0.025,0.034) 

(0.012,0.013,0.040,0.005) 

(1,1,1,1), (3.35,1,1,1) 

(1,1,1,1), (3.35,1,1,1) 

(1,1,1,1), (3.35,1,1,1) 

(9,1,0.01,11), (10,200,0.0001,8) 
(9,1,0.01,11), (10,200,0.0001,8) 

(19,11,0.001,10), (12,210,0.00001, 7) 


By selecting the scaling and shape parameters according to the conditions prescribed in Theorem 
15.11 (see Tabled]), a rational cubic spline FIF (cf. Eq. (13.21) 1 lying above the prescribed polygonal 
is generated in Fig. Oja)- We construct rational cubic spline FIF in Fig. dlb) by changing the 
scaling vector a with respect to the parameters of Fig. dla). By taking the null scale vector, 
i.e., a = (0,0, 0,0), we retrieve the classical rational cubic spline plotted in Fig. die). 

With specific choices of parameters (see Tabled! satisfying conditions prescribed in Theorem l5.21 
we obtain constrained rational cubic FIFs displayed in Figs, dla)-(c) lying above the quadratic 
spline q defined as follows. 


q(x) 


' X 2 -73+20 
2 ’ 

(x-3) 2 -(x-3)+8 
. 2 ’ 

— (x —7) 2 +21(x —7)+60 
6 > 

87(x —10) 2 —375(x —10)+200 
, 50 


if 0 < x < 3 
if 3 < x < 7 
if 7 < x < 10 
if 10 < x < 15 


In contrast to the classical rational spline f = f°, the rational cubic spline FIF /“ has derivative 
(/“)' having nondifferentiability in a finite or dense subset of the interpolation interval. Further, 
the irregularity may be quantified in terms of box-counting dimension, which depends mainly 
on the scaling vector a. This may find potential applications in various nonlinear and nonequi¬ 
librium phenomena wherein smooth interpolant satisfying a restraint (for instance, positivity, 
or more generally lying above or below a prescribed curve) and possessing irregularity in the 
derivative of suitable order is required. 
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FIF. 




(b) Effect of change in a in (c) Effect of change in shape 
Fig. Ufa). parameters in Fig. [2Ja). 


Figure 2. Rational cubic spline FIFs lying above a quadratic spline q. 


Note: Parts of the results of this paper were presented by the third author at the International 
Conference on Applications of Fractals and Wavelets (ICAFW-2015) held at Amrita School of 
Engineering, Amrita Vishwa Vidyapeetham, Coimbatore, Tamilnadu, India. 
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